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We derive a formula for the entropy of two dimensional incompressible inviscid flow, 
by determining the volume of the space of vorticity distributions with fixed values for 
the moments Q k = J u)(x) k d 2 x. This space is approximated by a sequence of spaces 
of finite volume, by using a regularization of the system that is geometrically natural 
and connected with the theory of random matrices. In taking the limit we get a simple 
formula for the entropy of a vortex field. We predict vorticity distributions of maximum 
entropy with given mean vorticity and enstrophy; also we predict the cylindrically sym- 
metric vortex field with maximum entropy. This could be an approximate description of 
a hurricane. 

1. Introduction 

Two dimensional inviscid incompressible flow is completely determined by so- 
lutions to the Euler equations in the absence of random external forces. Detailed 
numerical as well as analogue simulations!!] show complex phenomena such as recom- 
bination of vortices. Still, the situation is much simpler than in three dimensions. 
It appears that energy flows from small scales to large scales so that over long time 
intervals the flow becomes smoother. This is the opposite of the behavior observed 
in three dimensional hydrodynamics. Thus two dimensional flow is not turbulent in 
the same sense as three dimensional flow. Nevertheless, there is sufficient complexity 
remaining even in two dimensions that a statistical approach is worthwile. 

The study of the statistical mechanics of inviscid flow dates back to Onsagerl 
with some recent revivalsoQ. The crucial input into any statistical theory is the 
formula for entropy. Previous studies have used various postulates for the entropy 
of a vortex field u(x) such as 

ui(x) log \oj(x)\d 2 x. (1) 
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We will determine the formula for entropy predicted by its microscopic (Boltzmann) 
definition as the logarithm of the volume of the phase space with a fixed value 
of the conserved quantities (more precisely, central observables — see below). The 
conserved quantities of the Euler equation are the moments Qk — J u>{x) d 2 x. The 
volume of the phase space is difficult to determine directly because it is infinite 
dimensional. 

We will determine a formula for entropy by 'regularizing' the system; i.e., ap- 
proximating it with a finite dimensional system. The volume of the phase space 
is then finite — indeed the answer is known in the literature on Random Matrix 
theory!. The entropy of the original system can then be determined by taking the 
limit as the dimension goes to infinity. We establish this way that the formula for 
entropy is 

-LJlog\cj(x)-u(y)\d 2 x<?y (2) 

where A is the area of the region within which the flow is contained. This is 
quite different from the postulates used in earlier analysestrtJa. We then use this 
formula to predict the maximum entropy configuration with given mean vorticity 
and enstrophy* — it is the Wigner semi-circular distribution. We predict the vorticity 
distribution of an axi-symmetric vortex of maximum entropy. 

PolyakovQ has presented another theory of two dimensional turbulence, based 
on the group of conformal transformations. Our considerations are in many ways 
orthogonal, being based on the group of area preserving transformations rather than 
conformal transformations (which preserve instead the angles). 

In Section 2 we present an overview of two dimensional hydrodynamics. The 
Euler equations are formulated as a hamiltonian dynamical system, in analogy to 
the rigid body equations also due to Euler. 

In Section 3 we discuss the regularization procedure. We use ideas that have their 
origin in Quantum field theory, but have already been used in the hydrodynamics 
literature! 

In Section 4 we reformulate the regularized system in terms of hermitean matri- 
ces which allows us to use the ideas from random matrix theory to derive a formula 
for entropy in Section 5. In Section 6 we derive the distribution function for vortic- 
ity that maximizes entropy for given mean vorticity and enstrophy. And finally we 
predict the axi-symetric vortex field of maximum entropy in Section 7. 

2. Two-Dimensional Hydrodynamics as a Hamiltonian System 

Incompressible hydrodynamics is described by the Navier-Stokes equation: 

^ = -Vp + i?- 1 Au + f (3) 

along with the constraint V ■ u = 0. Here, u,p, R, f are the velocity field, the 
pressure, the Reynolds number, and the external force field respectively. We are 

'Enstrophy is the second moment of vorticity: Q2 = J ui 2 (x)d 2 x. 
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using units in which the mass density of the fluid is equal to one. Also, the total 
derivative is defined by 

Da da _ 

m-m + ^ a - (4) 

Pressure can be eliminated by taking the curl of this equation. This gives 
Du> 

— =w Vu + R~ 1 Au> + V x f, (5) 

where u) = V x u is the vorticity. 

Note that for incompressible flow, vorticity determines velocity: 

u a (x) = e abc d b J iv c (y)G(x,y)d 3 y, (6) 

or 

u a [x) = e ab d b / uj(y)G(x,y)d 2 y, (7) 



in two dimensions. Here G(x,y) is the Green's function of the Laplace operator 
A x G(x, y) = —S(x — y) with the appropriate boundary conditions on vorticity. For 
example, if u> — > as \x\ — ► oo, G(x,y) = — iv \ x-y \ ' n three dimensions and 

G(x, J/) = — 27r l°S I 2 - ~ 3/ 1 m t w0 dimensions. 

If the flow is two dimensional (i.e., u is independent of the z coordinate), vorticity 
is orthogonal to the ccy-plane so that u> ■ Vu = 0. Moreover vorticity can be 
thought of as a pseudo-scalar. Thus, for two dimensional incompressible inviscid 
(i.e., Br 1 — * 0) flow without external forces, vorticity is conserved along flow lines? 
with 

-p- =0, lu = e ab d a u b . (8) 

By expressing velocity in terms of vorticity we get an integro-diffcrcntial equa- 
tion: 

— = d a uje ab d b / G(x,y)uj(y)d 2 y. (9) 

These are analogous to another system studied by Euler: the equations for the 
rigid body given by 

— T~ — e abcL b A c dLd, (10) 

at 

where A ab is the inverse of the moment of inertia. The integral operator on the rhs is 
analogous to the inverse of moment of inertia. Thus vorticity corresponds to angular 
momentum, the stream function to angular velocity and the Laplace operator to 
the moment of inertia. The deep reason for this analogy! is that both systems 
describe equations for geodesies on a group manifold. In the case of the rigid body, 

"For three dimensional flow, vorticity is a vector field. Inviscid flow without external forces satisfies 
+ CxlIjJ = 0, where the Lie derivative of vector fields is C is given by Cuuj = u ■ Vu — ui ■ \7u. 
Thus, vorticiy is still conserved by the flow. 
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it is the rotation group with the metric given by the moment of inertia tensor; in 
the case of hydrodynamics, it is the group of area preserving diffeomorphisms with 
the L 2 -metric. In more physical terms, this analogy can be seen in the hamiltonian 
formalism. 

Recall that the space of functions on the plane is a Lie algebra, with the Lie 
bracket given by: 

[fi,f2]=e ab d a f 1 d b f 2 . (11) 

Every function / : R 2 — > R corresponds to an observable of two dimensional hy- 
drodynamics, uif = J f(x)uj(x)d 2 x. The above Lie bracket of functions suggests a 
natural postulate of Poisson bracket for vorticity: 

{uj(x), uj(y)} = e ab d b cu(x)d a S(x - y); (12) 

so that, {ojf 17 ujf 2 } = Ldy 1 j 2 ]. The natural postulate for the hamiltonian is the total 
energy of the fluid, 

H=\j u 2 d 2 x = X - J G(x, y)u{x)uj{y)d 2 xd 2 y. (13) 

A straightforward calculation shows that this hamiltonian with the above Pois- 
son bracket for vorticity indeed gives the equation of motion of two dimensional 
hydrodynamics. 

This is analogous to the hamiltonian formalism of rigid body mechanics. The 
Poisson bracket of the angular momentum (as measured in the co-moving coordi- 
nates) arises from the Lie algebra of infinitesimal rotations, 

{L a ,L b } = e abc L c . (14) 

The hamiltonian is the rotational kinetic energy 

H = l -A ab L a L b , (15) 

where A is the inverse of the moment of inertia matrix. Thus, vorticity is analogous 
to angular momentum; the Lie algebra of functions on the plane analogous to the 
Lie algebra of rotations; and the Green's function of the Laplace operator analogous 
to the inverse of the moment of inertia matrix. 

The square of the angular momentum L 2 = L a L a is a central function (i.e., 
has zero Poisson bracket with all functions of angular momentum) . In particular it 
commutes with the hamiltonian and hence is conserved. The phase space of rigid 
body mechanics is the sphere on which 1? is constant: this is the 'symplectic leaf 
of the Poisson algebra. 

Analogously, the moments 



Q k = / uj"{x)d 2 x 1 k = 1,2,- ■■ (16) 
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are central functions (hence conserved quantities) in two dimensional hydrodynam- 
ics; the phase space of two-dimensional hydrodynamics is the set of all vorticity 
fields with a given set of values of these Qk- 

The information in these moments can be packaged into the vorticity distribution 
function p(X), 



Geometrically, p(X) is the perimeter of the contour curve of vorticity where u){x) = 
A. This can be written explicitly as 



The shape of these curves changes with time but not their perimeter. 

In addition to these central functions, there may be additional conserved quan- 
tities such as momentum and angular momentum: 



But the presence of a boundary can break translational or rotation invariance, 
violating these conservation laws. 

In spite of the presence of these infinite number of conservation laws, two dimen- 
sional hydrodynamics is far from being integrable: the phase space on which the Qk 
are constant is still infinite-dimensional. The crucial step in a statistical approach 
is the correct identification of entropy. This is ultimately defined by the canonical 
structure of the phase space: entropy of a macroscopic configuration is the log of 
the volume of the phase space corresponding to it. We would for example like to 
determine the volume of the phase space with a given set of values of moments 
Qk- However, this is some infinite dimensional subspace of the space of all vorticity 
distributions. We need to find a way to approximate the phase space by a finite 
dimensional space: a 'regularization' or 'discretization' of the system. The entropy 
can then be determined within this discretized version of the theory and eventually 
a limit as the number of degrees of freedom goes to infinity can be taken. 

We now describe an elegant discretization! of two dimensional hydrodynamics 
that preserves its symmetries (Lie algebra structure) and conservation laws. We 
will see that the problem of determining the entropy of a two dimensional flow can 
then be solved using ideas from Random matrix theory due to Wigner. 

3. Regularization 

It will be convenient to assume that the flow satisfies periodic boundary con- 
ditions in the plane. Our final result for entropy will be the same even if other 
boundary conditions are assumed, but the intermediate formulas seem simplest 
with periodic boundary conditions: 




(17) 




(18) 




(19) 



(20) 



6 Two Dimensional Turbulence 



We can then Fourier analyze vorticity: 



Li 



uj(x)e 



' d 2 x, 



(21) 



uj{x) 



1 



L1L2 



E 



UJt, 



g Lj g i 2 



mi ,m2 — — 00 

The hamiltonian and Poisson bracket are, in terms of these variables, 

2tt 



H = (L\L 2 f V* -\\uj m \ 2 , {ci) TO ,cJ n } = - 

z — ' TCI 



m#(0,0) 

The equations of motion become 



-t a bm a nbU} m+n . 



(22) 



(23) 



dWp V"^ ab 

~df = E £ ™- 

m+n— p 



H6 



1 1 



W m U) n - 



(24) 



We are still dealing with a system with an infinite number of degrees of freedom. 
At first glance, ignoring all except the low momentum modes (i.e., keeping only 
I Til |, Imal < N for some N) looks like a reasonable 'regularization' of the problem: 
we would lose only information at very small length scales. However, such a naive 
truncation that simply ignores large mi, ma modes would not be consistent: for 
example, the Poisson bracket of two low momentum modes could be a high mo- 
mentum mode. We must modify the Poisson brackets so that the Poisson brackets 
of the modes we keep form a closed Lie algebra; as N — > 00 this modification must 
dissappear and we must recover the original Lie algebra. Moreover the hamiltonian 
must be modified so that the equations don't mix high and low momentum modes. 

A formalism along these lines is known in the literatureS. The essential idea is 
to modify the coefficients in the Poisson brackets and hamiltonian so that they are 
periodic with period N; this way the algebra would be closed. There is a choice 
that satisfies the Jacobi identity, thus preserving the Lie algebra structure: 



1 



{u> m , uj n } = 7j sin[0(min 2 - m 2 n.i)]w m+ „ mod jv, 6 = 



2?r 
TV' 



(25) 



We keep only a finite number of Fourier modes (we can assume that N is odd 
N = 2v + l), 

mi, m2, ni, ^2 = — v, ■ ■ ■ 0, 1, • • • v. (26) 
To be fair, the hamiltonian is truncated as well: 



1 



1 



H 2 E A ( m )l 
mi, ni2 = —v 
A(m) ^ 



(27) 
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It is natural to deform the eigenvalues of the Laplacian to A(m) = { ^ sin [^mj } + 
{^-j sin [^012] } 2 . This would preserve the periodicity of the eigenvalue modulo N. 

It is clear that as N — > 00, 9 — ► and the structure constants of the Lie 
algebra tend to the original ones; so does the hamiltonian. The Lie algebra we 
obtain this way is in fact nothing but that of the Unitary group U(N), the Lie 
algebra of hermitean matrices. Indeed the discretization procedure above has a 
natural interpretation in terms of non-commutative geometry; the Lie algebra is 
the algebra of derivations of the algebra of functions on a non-commutative torus. 
The constants A(m) in the hamiltonian above are the eigenvalues of the Laplace 
operator. However, we will not need this interpretation in what follows. Non- 
commutative geometry plays a deeper role in the three dimensional version of this 
theory which one of us is still developing. 

4. Matrix Formulation 

It will be convenient to make a linear change of variables that will make the 
connection with hermitean matrices more explicit. Define N x N unitary matrices 
Ui , U2 satisfying 

C/iC/ 2 = e ie U 2 U!. (28) 

Then, U{m) = e-i^^U^Up satisfies W{m) = U{-m). To be specific we 
choose, 




... e i{N-\)0 1 



\ /0 ••• 1 

1 ••• 



(29) 



\0 ••• 1 0, 



(These matrices can be viewed as the coordinates on a non-commutative torus.) 
Now the Fourier coefficients of vorticity can be packaged in the hermitean matrix 

Co = u m U(m). (30) 

m 

(Recall that lo(x) being real implies that uj* n — u)- m . This along with U^(m) — 
U(—m) implies that = Co.) Now we can verify that 

{uj ab , Cj cd \ = i[5 bc uj ad - 6d a u; c b\. (31) 

These are the well-known commutation relations of the Lie algebra of U(N) in 
the Weyl basis; this establishes the identification of our truncated Lie algebra. 
Moreover, we can see that the quantities 

Q fe = tiLu k , k = 0,l,---N- 1 (32) 

are central in this algebra: they commute with any function of Co. In the limit as 
N — > oo, 9 — > 0, these tend to the central functions Qk we started with. 
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The hamiltonian becomes, 

H = -Wab&cdGabcd (33) 

for some tensor G a b c d whose explicit form we do not need. The equation of motion 
are 

^=mM (34) 

where we have Q a b = G a bcd^>cd- Thus time evolution is a sequence of unitary 
transformation so that the traces of powers of u> are unchanged. 



5. Formula for Entropy 

Now it is clear that the information contained in the moments Q\ — tr <2> k = 
Yjj=i 1S simply the spectrum {Ai, ■ • • , Ajv} of the hermitean matrix ui. The phase 
space of the system is a vector space of finite dimension N 2 . The submanifold of 
matrices with a fixed spectrum {Ai, • • • A at} is compact and has dimension N(N — 1): 
it is called the 'flag manifold' in algebraic geometry. There is a unique (up to a 
multiplicative constant) volume form on this submanifold invariant under the action 
of the unitary group. (This is the volume form induced by the symplectic structure 
associated to the Poisson brackets above.) The volume of microstates with a fixed 
value of moments Qk are given by integrating thia volume form. 

The volume of this submanifold is well knownO: 

V N = C N I] l x k-M 2 - (35) 

l<fc<i<JV 

Here CV is a constant whose value we do not need. Thus the log of the volume of 
this manifold becomes 

\ogV N = 2 log|A fe - A^+logCjv. (36) 

l<fc<(<JV 

In terms of the density distribution of eigenvalues 

1 N 

/c=i 

this becomes 

logV N =N 2 V J PN {X) PN {X')log\X- X'\dXdX' + \ogC N . (38) 
Here V J denotes the principal value integral: 

V [ /(A, X')dXdX' = lim / f(X, X')dXdX' . (39) 

J e~>0+ J\ x _x,\ >e 



Using Random Matrix Theory 9 



The advantage of this point of view is that it survives the limit N ~ * oo, except for 
an undetermined additive constant. 
Now, in the limit N — > oo, 

p N (X)X k dX^ jJp(X)X k dX, (40) 

where p(X) — J S(uj(x) — X)d 2 x as we noted earlier. Also A = J d 2 x is the total 
area of the region of the flow. Thus we see that limjv^oo Pw(A) = /5(A), where 
p(X) = \p{X) is the normalized distribution of vorticity: J p(X)dX = 1. We can 
now get the limiting form of the log of the volume of the submanifold with a given 
set of moments: 

XH ^lim^logV^^ J p(X)p(X') log |A - X'\dXd\' . (41) 

(We drop the undetermined additive constant.) We can rewrite this directly in 
terms of vorticity: 

d 2 xd 2 y 



XH = V j log \uj{x) - U (y)\—p±. (42) 

This quantity should be viewed as the entropy of a vorticity distribution. 

This formula for entropy is one of our main results. Our careful derivation by 
truncating the Lie algebra gives an entirely different result from those postulated 
in earlier analysesoiu, namely, 

Lu{x)log\uj(x)\d 2 x. (43) 



6. Semi-circular distribution 

As an application of this new notion of entropy we describe here the vortex 
distribution of maximum entropy with given total vorticity Q\ and enstrophy Qi. 
We need to maximize 

X(P) = V [ p(X)p(X') log | A - X'\dXdX' (44) 



keeping the quantities J p(X)X k dX for k = 0,1,2 fixed. Introducing Lagrange mul- 
tipliers Co, ci, C2 for the constrained variational problem, we get 

2V J p(A') log |A - X'\dX' = c + ciX + c 2 A 2 . (45) 

Differentating this to eliminate Co gives the singular integral equation 

2V j Y^Ji dx ' = ci + 2c 2 A. (46) 
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The equation is well known in the theory of random matrices; the solution is the 
'semi-circular' distribution of Wignem 

PW = (|A - Qx\ < 2a) V [4a 2 - (A - Qi) 2 ] . (47) 




lambda 



Fig. 1. Vorticity distribution function of maximum entropy with Qi = 5.0, a = 1.0. 



Here Qi — \ f ui(x)d 2 x is the mean value of vorticity. The constant a is the 
standard deviation of the vorticity distribution and is determined by the enstrophy 
per unit area Q2 = \ J u>(x) 2 d 2 x: 



a 2 . 



(48) 



7. Vorticity Distribution in a Cylindrically Symmetric Fluid 

As an applicaton of our ideas we determine the maximum entropy configuration 
of an axially symmetric fluid trapped between two concentric circles of radii a± < 02- 
The vorticity is a function only of the distance from the origin in the plane. 

With axial symmetry, any function uj(r) of the radial distance alone will be a 
static solution of the Euler equations. We will seek the answer to the following 
question: among all cylindrically symmetric functions with given mean vorticity 
and enstropy, which one has the largest entropy? One should expect that a system 
subject to many perturbations will eventually settle down to this distribution. 

Now, the vortex contours of an axisymmetric vorticity distribution are concentric 
circles. We should expect that the vorticity is monotonic with distance in the stable 
configuration. The vortex distribution is given by 



P(A) 



S(p(r) - \)d z x = 2tt / 5{uj(r) - X)rdr 



(49) 
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Let us introduce the variable 



Rather than think of the vorticity w as a 



function of r (or u), let us try to determine its inverse function; that is, it as a 
function of u. Then we see that 



p(X) = 7T 



du 
dX 



(50) 



From the last section we know p(X) = Ap(X) to be the semi-circular distribution. We 
can solve the above differential equation to determine the corresponding vorticity 
profile in parametric form: 



uj = 2cr sin ( 



IT 



i sin (20) 



with the parameter <j> taking the range of values 



2 ~ y 2 



(51) 



(52) 



There are two possible solutions with equal entropy: one with vorticity increasing 
with distance and the other with it decreasing. In most situations we would ex- 
pect the vorticity to decrease with radial distance. We give a sample plot of this 
distribution below. 




10 12 14 



radius 



Fig. 2. Vorticity field of maximum entropy within radii a\ = 1.0, <J2 = 15, mean vorticity Qi = 5.0, 
and variance ct 2 =1.0. 



This could describe a vortex such as a hurricane, Jupiter's red spot or a tornado, 
to a reasonable approximation. A comparison with experimental measurements 
would be interesting. 
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